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ABSTRACT 

Three types of turbulence models which account for rotational effects in noniner- 
tial frames of reference are evaluated for the case of incompressible, fully developed 
rotating turbulent channel flow. The different types of models are a Coriolis-modified 
eddy-viscosity model, a realizable algebraic stress model, and an algebraic stress 
model which accounts for dissipation rate anisotropies. A direct numerical simu- 
lation of a rotating channel flow is used for the turbulent model validation. This 
simulation differs from previous studies in that significantly higher rotation numbers 
are investigated. Flows at these higher rotation numbers are characterized by a re- 
laminarization on the cyclonic or suction side of the channel, and a linear velocity 
profile on the anticyclonic or pressure side of the channel. The predictive performance 
of the three types of models are examined in detail, and formulation deficiencies are 
identified which cause poor predictive performance for some of the models. Criteria 
are identified which allow for accurate prediction of such flows by algebraic stress 
models and their corresponding Reynolds stress formulations. 

I INTRODUCTION 

Turbulent flows in noninertial reference frames are of considerable interest in a 
variety of industrial applications. However, the success of a computational analysis 
of such flow phenomena relies heavily on the choice of turbulence model. It is well- 
known that, without some modifications, conventional isotropic eddy-viscosity models 
fail to predict the effect of noninertial forces on turbulence, whereas second-moment 
closures, for example, can account for noninertial effects in a systematic way. 

An explicit algebraic stress model (EASM) has been developed by Gatski and 
Speziale (1993) which is a nonlinear extension to an isotropic eddy-viscosity two- 
equation model. This EASM allows for the inclusion of stress anisotropies and rota- 
tion rate effects and is developed on a rigorous mathematical basis, while still keeping 
the associated computational cost comparable to that of two-equation models. Pre- 
vious results have shown that this approach is a viable approximation to the full 
Reynolds stress closure for a variety of engineering flows. An extension to this EASM 
which accounts for dissipation rate anisotropies in a systematic way and is applicable 
to wall-bounded flows, has been recently evaluated (Xu and Speziale 1996), and is 



based on an analysis of the exact transport equation for the dissipation rate tensor 
(Speziale and Gatski 1997). This composite algebraic stress model (CASM) is ex- 
tended here to include noninertial effects in the determination of both the stress and 
dissipation rate anisotropies. The composite model differs from previously proposed 
algebraic stress models in two ways: dissipation rate anisotropies are accounted for 
in the constitutive relation and the coefficient of the production term in the trans- 
port equation for the scalar dissipation rate is now sensitized to the mean strain and 
rotation rate tensors. 

The objective of this study is to systematically explore the predictive capabilities 
of explicit algebraic stress models in general and, in particular, the composite model 
for the case of a fully developed turbulent channel flow with strong spanwise rotation. 
This effort is an extension of previous work in that (1) the recent DNS obtained by 
Lamballais et ai (1996) has yielded new data for rotation numbers that are higher 
than those for the previous DNS study of Kristoffersen and Andersson (1993) (as well 
as the previous experimental study of Johnston et ai.1972) and (2) the equations are 
integrated up to the wall, whereas most of the previous computations were limited to 
wall-function boundary conditions that are not suitable for (strongly) rotating flows 
where regions of relaminarization can occur. 


II TURBULENCE MODELS 

The incompressible, fully developed rotating channel flow is a unidirectional flow 
that results in simplified expressions for the mean flow in the noninertial frame. In 
this frame, the mean strain rate and rotation rate tensors 
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where S = S(y) and y is the distance measured from the (bottom) wall of the channel, 
as shown in Figure 1. 

The common feature between linear and nonlinear eddy viscosity models is both 
require the solution of only two transport equations; an equation for the turbulent 
kinetic energy h 
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Figure 1. Schematic of fully developed turbulent channel flow in a rotating frame 


V — — 2tuS is the turbulent production, r r2 is the Reynolds shear stress, f 2 and / M 
are wall damping functions, is a closure constant, and v is the kinematic viscosity. 


Three turbulence models will be evaluated by using (3) through (5): a Coriolis- 
modified eddy- viscosity model (EVM) proposed by Howard et al. (1980); , which is 
based on the Launder and Sharma (1974) model; a realizable algebraic stress model 
(ASM) (Shih et al. 1995); and a CASM that accounts for dissipation rate anisotropies 
(Xu and Speziale 1996). These three models are distinct in the way rotational effects 
are incorporated into their formulations. 

In the EVM, the noninertial modification is phenomenological and based on an 
analogy with curved boundary layers. This type of modification was used recently by 
Pettersson et al. (1996) to rotating channel flow but at much lower rotation numbers. 
The realizable ASM model of Shih et al. (1995) includes the effect of mean rotation 
and has been tested on rotating homogeneous shear flow. The CASM is an algebraic 
stress model which is directly extracted from a full Reynolds stress closure as well as 
a transport equation for the tensor dissipation rate £,-j. Thus, the model includes the 
effects of an anisotropic dissipation rate in a algebraic stress framework. 

The Coriolis-modified eddy-viscosity model uses the Boussinesq-type relation for 
the Reynolds stresses: 
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where C; = C„ = 0.09, Q = 1.44, C e2 - 1.92, a K = 1, cr e = 1.3, and D is the 
rotation rate of the reference frame. 

The Reynolds stress algebraic equation model (Shih et a 1. 1995) is given by 
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and 
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t/ r is the friction velocity, — 1.0, = 0.09, C*i ~ 1-44, C* 2 = 1.92, a/c = 1, and 

(7 5 = 1.3. 

The composite algebraic stress model that accounts for dissipation rate anisotropies 
(Xu and Speziale 1996) is a extension of the EASM of Gatski and Speziale (1993) in 
inertial frames, in the noninertial case, the process is complicated by the fact that 
the system rotation enters differently into the anisotropic dissipation rate and the 
algebraic stress relations. 

First, noninertial effects are introduced into the explicit algebraic anisotropic dis- 
sipation rate equation developed by Speziale and Gatski (1997) as 
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( V, = 5.8. and = 0.6. 

In the absence of rotation, this explicit algebraic anisotropic dissipation rate model 
could be directly inserted into the explicit algebraic stress model and the resulting 
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composite explicit algebraic equation could be solved in conjunction with the tur- 
bulent kinetic energy and the dissipation rate equations. In the noninertial frame 
this direct substitution is not possible, and the problem can be quickly identified 
by examining the implicit relation for the algebraic stress model that is used as the 
starting point for the tensor polynomial expansions associated with the explicit rep- 
resentations. The implicit relation that accounts for dissipation rate anisotropies can 
be written as 
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Ci = 3.4 + 1.8'P/c, C -2 = 0.36, C3 = 1.25, and C 4 = 0.40. These constant closure 
coefficients are obtained from the SSG pressure-strain correlation model (Speziale et 
a 1. 1991). A comparison of (17) and (22) clearly shows that in general the rotation rate 
tensors in the noninertial frame that are extracted from the algebraic dissipation rate 
model and the algebraic stress model are not the same. This difference precludes a 
simple combination of terms as suggested by (21) and shown to be possible in inertial 
frames (Xu and Speziale 1996), where integrity bases were used to get a composite 
explicit representation for the Reynolds stresses. 

As an alternative to this approach, consider (21) rewritten as the matrix system 
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The system in (24) can be inverted analytically to obtain explicit expressions for 
the Reynolds stress anisotropies i> tJ (and tv,). Once again, these stress relations are 
coupled with (3) and (4) for the turbulent kinetic energy and the turbulent dissipation 
rate, with 
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f n = 1.0, (7 M = 0.094, C* 2 = 1.83, a k = 1, and a e = 1.3. The additional (noncon- 
stant) term in the expression for C* x represents a production ratio of the turbulent 
dissipation rate {d^S) and turbulent kinetic energy (6 12 5'). 


Ill RESULTS 


The turbulence models presented in the last section are coupled with the stream- 
wise momentum equation 


du 2 dS dr l2 
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where «, = ((/, 0,0); G is the (constant) effective pressure gradient, which includes 
the centrifugal force term; and the bulk Reynolds number Reb = hUb/v, with l ! b 
as the bulk velocity. Equation (30), coupled with the transport equations for the 
turbulent kinetic energy (3) and the turbulent dissipation rate (4), is integrated to 
steady state by a one-dimensional second-order finite-difference scheme. This simple 
one-dimensional spatial problem allows for solutions with arbitrarily high numerical 
accuracy by using a sufficient number of points. Here, 200 points were typically used, 
with highly stretched meshes near the solid walls. 

Because the different rotation regimes considered by the DNS (Lamballais et al. 
1996) were obtained at the same bulk Reynolds number, the pressure gradient in the 
numerical code was adjusted in order to have f?e& = 2500 at convergence. In the 
following, results for three different rotation numbers, defined by 

Ro = 2\n\h/U h (31) 


are shown (Ro = 0, Ro = 0.5, and Ro = 1.5). In the earlier DNS study of Kristof- 
fersen and Andersson (1993), the maximum rotation number studied was 0.5. Thus, 

the DNS data used in this study significantly increases the validation range for the 
turbulence closure models. 

The profiles of the mean velocity and turbulent kinetic energy are shown in Figures 
2 and 3. respectively. The turbulent kinetic energy, as well as the turbulent stresses 
to be presented, are scaled by an average friction velocity u T , which is the half-sum of 
the friction velocities on both walls. The DNS results show the characteristic linear 
region of slope 2D in the mean velocity, which leads to a mean absolute vorticity 
2 (S' — Q) that is close to zero. These results also show that, the turbulent kinetic 
energy is higher on the anticyclonic or pressure side (y = 0) than on the cyclonic 
or suction side (y = 2), where relaminarization occurs. In Figure 2, for the mean 
velocity the composite model is able to reproduce all features of the flow for the 
three rotation numbers considered, including the linear portion of the profile and the 
relaminarization on the cyclonic side, characterized by a parabolic velocity profile. 
For the turbulent kinetic energy shown in Figure 3, the asymmetry of the profile 
and the higher turbulence intensity on the anticyclonic side of the channel are clearly 
visible and are consistent with the DNS results. In the case for which Ro = 0, the 
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peaks in turbulent kinetic energy near the wall are not well predicted because in the 
CASM no ffj damping function has been introduced. 

The Coriolis-modified model of Howard et a/. (1980) gives reasonable predictions 
for the mean velocity at the different rotation numbers although not with the same 
degree of accuracy as the CASM. The algebraic model of Shill et al. (1995) does 
not correctly predict the mean velocity at these higher rotation numbers since it 
fails to predict the linear profile on the anticyclonic side and the relaminarization 
on the cyclonic side. For the turbulent kinetic energy, both the Howard and Sliih 
models misrepresent the behavior of the kinetic energy in the higher rotation rate 
cases. At Ro = 0.5 and Ro = 1.5, the Coriolis-modified EVM yields results that 
effectively damp out the turbulent kinetic energy in the relaminarizing portion of the 
flow, whereas the Shih algebraic stress model is somewhat insensitive to the effects 
of rotation. 



Figure 2. Mean velocity profiles for (a) Ro = 0. (b) Ro = 0.5, and (c) Ro = 1.5. 

A — A DNS (Lamballais et al. 1996); . CASM; , ASM (Shih et al. 1995); . 

Coriolis-modified EVM. 





Figure 3. Turbulent kinetic energy profiles for (a) Ro = 0, (b) Ro = 0.5, and (c) 

Ro = 1.5. A — A DNS (Lamballais et al 1996); , CASM; , ASM (Shih et al 

1995); , Coriolis-modified EVM. 


With the success of the CASM, further investigation of the flow dynamics is worth- 
while by examining the total shear stress. At steady state, (30) can be integrated 




(32) 


with respect to y and expressed in wall units as 


*12 + 


Z 9+ _ 0,2 

Re r 



where u t0 = u T | y=0 , u T i = u T | y= 2 , and Rc T — hu r jv. As this equation shows, the 
total shear must vary linearly across the channel for all rotation rates. The partition 

of the total shear stress between the turbulent and viscous 2 S + / Re r stresses is 
illustrated in Figure 4 for the composite model. The region of neutral stability, where 
the velocity profile is linear, is characterized by a turbulent shear stress that varies 
linearly, and a viscous shear stress that remains constant. 



y/h y/h y/h 


Figure 4 . Partition of the total shear between turbulent and viscous shear stresses: 

computations by ('ASM for (a) Ro = 0, (b) Ro = 0 . 5 , and (c) Ro = 1.5. , 

+ 2 S + /Re T ; , -r^; , 2 S' + /Rt T . 

In addition to the turbulent shear stress, the normal Reynolds stresses are signifi- 
cantly affected by the rotation. Figure 5 shows the attenuation of the normal stress 
components on the relaminarized side of the channel at both nonzero rotation num- 
bers. On the turbulent side, the streamwise component Tn is attenuated relative to 
the Ro = 0 case, and both the T22 and 733 components are enhanced relative to the 
Ro — 0 case. The most significant effect is on the r 2 2 component, but even at the high 
rotation case the 733 component also exceeds the T\\ component on the anticyclonic 
side. This result is consistent with the DNS results and can be explained by the fact 
that the ( 4 ASM is derived from a Reynolds stress model and will, therefore, inherit 
the right sensitivity of the production terms for the individual normal stresses to the 
rotation. On the other hand, the stress anisotropies predicted by the model of Shill 
et al. ( 1995 ) in ( 11 ) for unidirectional shear, where b n — —622 = 2 c 2 (*Sr) 2 (l — -) and 
633 — 0, show that this coupling of the b u and 6 22 forces the incorrect prediction of 
isotropic turbulence when 0 = , 9 , and which would preclude the correct prediction of 
the normal stresses as displayed in Figure 5 . 

Both the presence of the linear velocity profile and the relaminarization process 
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Figure 5. Turbulent normal stress profiles for (a) Ro = 0. (b) Ro = 0.5, and (c) 

Ro — 1.5. DNS (Lamballais et a 1. 1996) data: A — A, Tn; □ □, T 22 ; o o, 733 . 

CASM: , Tu - - - T 22 j , T33- 


can be explained. The CASM is derivable from the Reynolds stress model 
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where P tJ (and TV = 'Dj j/2) represent the effects of turbulent and viscous diffusion; 
oi, a 2 , O 3 , and g are given in (23). The implicit form of the algebraic stress model 
is obtained by setting the left side of (33) to zero and using the appropriate closure 
model for the dissipation rate anisotropies (such as the one shown in (16), which is 
used here). Note that the equilibrium hypothesis that underlies the algebraic stress 
models is exactly satisfied here (fq- = I\ = 0 ), and the only approximation that is 
made is TV, = jfPh ■ 

An examination of (33) shows that the anisotropy component b i2 must satisfy the 
following equation: 
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where 7Z is given in (27), and 
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As could be expected, this expression for b r2 has the same functional form as the 

EASM: 


b v2 — — Q]C*5r, C* 
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(36) 

Thus, in both the composite model and the EASM the production-to-dissipation rate 
ratio is always given by 


— = - 4 b u Sr = 4 q 1 C*( 5 't) 


( 37 ( 


9 





showing that — is a function of both 0/5’ and St. By interchanging this dependency, 
the behavior of 0/5 as a function of j and 5r, 
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can be studied. For comparative purposes, a corresponding fl/S relationship can also 
be obtained from an expression equivalent to (37) for the Shih et a 1. (1995) model; 
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Figure 6 shows the evolution of D/5 across the channel for the different rotation 
regimes and the three models considered. For the two rotating cases, several features 

are apparent. Starting from values near zero on the anticyclonic side (because 5 
is high near the wall), the DNS results clearly show a plateau at fl/S = 1, then 
quickly grow and change sign at the location of the maximum velocity ( S = 0). On 
the cyclonic side, D/5(< 0) then approaches zero with a y~ l behavior ( 5 ~ — y 
in the relaminarized region). The figure shows that the CASM closely follows the 
DNS results and accurately predicts the location of maximum velocity. With (38), 



Figure 6. Scaled rotation rate variation across channel for (a) Ro = 0, (b) Ro = 0.5, 

and (c) Ro = 1.5. A — A DNS (Lamballais et al. 1996); , CASM; , ASM 

(Shih et ai. 1995); , Coriolis- modified EVM. 

we can explain why the CASM model is able to predict such features, and we shall 
see that the reason is closely related to the expression for C*, which must show the 
correct dependency on 0/5 and Sr, By their nature, the CASM, and more generally, 
all of the algebraic stress models that are consistently derived from Reynolds stress 
models inherit the correct behavior for C*, while algebraic stress models that provide 
a C* expression based solely on constraints such as realizability, calibrations, arid 
phenomenological arguments may not have the correct behavior and will fail to predict 
the neutral stability region arid the relaminarized zone that is observed in the rotating 
channel. 

For larger values of St (i.e., > 3), the last term under the root in (38) is negligible 
compared with the other terms, and f 1/S becomes a function of j only, and takes 
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values on the two limit branches that correspond to the sign of S(= ±|S'|). Now, it 
is also easy to verify that these two values of f l/S rapidly become independent of an 
increasing j, and asymptote to either (Q/S)+ ( S > 0 ) or (fl/S)~ (S < 0). Equation 
(38), therefore, shows that for a wide range of values of Sr and 7 the value of Q/S 
becomes effectively independent of these parameters, and takes values close to (fi/S)*. 
These two limiting values only depend on the values of the model coefficients a/s, 
and have the following values for the SSG pressure-strain model: (C2/£) + = 0.992 
and (tt/S)~ = —0.103. In the channel away from the walls, the diffusion of I\ may 
be expected to be small, and we should have 7 ^ 1 and *SY>3. In this case, the 
scaled rotation rate will have values on the limit branches that will be close to the 
limit values (Q/S)*. 

Figure 7 illustrates this phenomenon by showing the evolution of the scaled rotation 
rate correlated with the variation of 7 across the channel. In Figure 7(a), the two 



Figure 7. Scaled rotation rate variation versus 7 in the channel for (a) CASM, (b) 

ASM (Shih et al. 1995), and (c) Coriolis-modified EVM. , Ro — .5; . Ro — 1.5; 

, limit branches (shown for (a) only). 


limit branches given by (38) with values of Sr > 3 are represented by dashed lines. 
At 7 = 0 and y = 0 , the curves that correspond to the CASM simulation at the 
two rotation regimes first move in a region where (Q/S)~ < Q/S < ( U / 8 ’ ) + because 
the values of Sr that are given by the model are also very small. However, as 7 
rapidly increases with movement away from the near-wall region, St increases also, 
and the points collapse on the (positive) limit branch. Until very near the wall at 
y — 2 , the value of St stays at values sufficiently high to force the points to stay 
on the limit branches. After having attained values close to (Q/S)+ — 0.992 for the 
major portion of the channel on the anticyclonic side, the maximum velocity point is 
reached. 7 becomes very small (as S — > 0 ) and the negative limit branch is followed 
after the maximum velocity. For most of the values il/S < (il/S )~ , the corresponding 
value of j on the negative branch is very small, and according to (3), we can expect 
dK/cH < 0, which leads to relaminarization in this region. 

This explanation for the occurrence of the linear profile and the relaminarization 
process is also valid for other Reynolds stress models (i.e. other pressure-strain cor- 
relation models). For example, the LRR model (Launder et al. 1975) yields limiting 
values for tl/S of (fl/S) + = 0.644 and (fi /S)~ = —0.131 for the anticyclonic and 
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cyclonic sides, respectively. However, one would expect that the slope of the linear 
velocity profile, when compared to DNS, would not he correct. Note also that these 
results are independent of the particular model for e that is used, because the effect 
enters only through Sr and 

The same analysis can be carried out on the Shih et a 1. (1995) model by using (39). 
With this model, the evolution of St always has an effect on Sl/S\ and no limiting 
behavior occurs. Therefore, the system is then not forced to reach the asymptote 
fl/S ~ 1, which precludes attainment of a linear velocity profile on the anticyclonic 
side (Figure 7 (b)). From ^ = 0 at j/ = 0, 7 quickly increases and the resulting 
curve is obtained from the balance of the model equations and (39). The fact that 
this model cannot reproduce the correct features for the rotating channel at high 
rotation numbers can be attributed to the lack of a mechanism in C* to render Sl/S 
independent of St and 

Finally, for the Coriolis-modified EVM, relation (37) cannot be used to find a 
relation between fi/5, and Sr. Instead, (4) can be examined at steady state in 
regions away from the walls where the damping functions and diffusive terms can 
be neglected. The dissipation rate equation then yields the simple relation that the 
production-to-dissipation rate ratio is C* 2 /C*i and 

n _ 1 

5 ” 2 

For sufficiently high values of Sr, the dependency on j is totally removed, and f l/S 
takes a value of 0 or 1. Other values of fl/S can only be reached when 7 - vanishes, 

as shown in Figure 7(c). Thus, the model then yields an abrupt and total damping 
of the turbulence. (See Figure 3.) 

IV CONCLUSIONS 

This study has shown that algebraic stress models consistently derived from Rey- 
nolds stress models inherit the correct dependency to rotation, and noninertial effects 
are automatically accounted for in a rigorous way. On the other hand, algebraic stress 
models that try to generalize the eddy- viscosity hypothesis in a phenomenological way 
are not necessarily directly extendible to non-inertial frames. As these results have 
shown, the nonlinear eddv-viscositv function C* must be constructed with the correct 
dependency on St and il/S. This study also demonstrated that the key features of 
the rotating channel flow were controlled by mechanisms only remotely linked to 
the dissipation rate equation, which would mean that phenomenological models that 
attempt to account for noninertial effects solely through modification of the source 
terms in the dissipation rate equation may not be properly accounting for essential 
flow physics. 
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